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O ; Abstract 

^ ' The interpretation of experimental results at RHIC and in the future also at LHC requires 

I ' very reliable and realistic models. Considerable effort has been devoted to the development 

Oh! of such models during the past decade, many of them being heavily used in order to analyze 

^ ' data. 

There are, however, serious inconsistencies in the above-mentioned approaches. In this paper, 
we will introduce a fully self-consistent formulation of the multiple-scattering scheme in the 
. framework of a Gribov-Regge type effective theory. 

5t , With the start of the RHIC program to investigate nucleus- nucleus collisions at very high energies, 

there is an increasing need of computational tools in order to provide a clear interpretation of the 
data. The situation is not satisfactory in the sense that there exists a nice theory (QCD) but 
we are not able to treat nuclear collisions strictly within this framework, and on the other hand 
there are simple models, which can be applied easily but which have no solid theoretical basis. A 
good compromise is provided by effective theories, which are not derived from first principles, but 
which are nevertheless self-consistent and calculable. A candidate seems to be the Gribov-Regge 
approach, and - being formally quite similar - the eikonalized parton model. Here, however, severe 
inconsistencies occur, as discussed in ||^. 

As a solution of the above-mentioned problems, we present a new approach which we call 
"Parton-based Gribov-Regge Theory": we have a consistent treatment for calculating (topo- 
logical) cross sections and particle production considering energy conservation in both cases; in 
addition, we introduce hard processes in a natural way. 

The basic guidehne of our approach is theoretical consistency. We cannot derive everything 
from first principles, but we use rigorously the language of field theory to make sure not to violate 



1 



basic laws of physics, which is easily done in more phenomenological treatments (see discussion 
above) . 

Let us first introduce some conventions. We denote elastic two body scattering amplitudes as 
T2^2 and inelastic amplitudes corresponding to the production of some final state X as T2-,x 
(see fig. ||). As a direct consequence of unitarity on may write the optical theorem 2ImT2^2 = 





Figure 1: An elastic scattering amplitude T2_»2 (left) and an inelastic amplitude T2—,x (right). 

'^^{T2->x){T2-,x)* ■ The right hand side of this equation may be literally presented as a "cut 
diagram", where the diagram on one side of the cut is {T2^x) and on the other side {T2^x)* ■, as 
shown in fig. ||. So the term "cut diagram" means nothing but the square of an inelastic amplitude. 







Figure 2: The expression y (T2_x)-(T2^x)* which may be represented as a "cut diagram". 

summed over all final states, which is equal to twice the imaginary part of the elastic amplitude. 

Before coming to nuclear collisions, we need to discuss somewhat the structure of the nucleon, 
which may be studied in deep inelastic scattering - so essentially the scattering of a virtual photon 
off a nucleon, see fig. 11 Let us assume a high virtuality photon. It is known that this photon 




Figure 3: Deep inelastic scattering: the general diagram (left) the "sea" contribution with a soft Pomeron at the 
lower end (middle) and the "valence" contribution (right). 



couples to a high virtuality quark, which is emitted from a parton with a smaller virtuality, the 
latter on again being emitted from a parton with lower virtuality, and so on. We have a sequence 
of partons with lower and lower virtualities, the closer one gets to the proton. At some stage, some 
"soft scale" scale Qg must be reached, beyond which perturbative calculations are no longer valid. 
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So we have some "unknown object" - indicated by "?" in fig. ^ - between the first parton and 
the nucleon. In order to proceed, one may estimate the squared mass of this "unknown object", 
and one obtains doing simple kinematics the value Ql/x, where x is the momentum fraction of 
the first parton relative to the nucleon. Therefore, sea quarks or gluons, having typically small 
X, lead to large mass objects - which we identify with soft Pomerons, whereas valence quarks 
lead to small mass objects, which we simply ignore. So we have two contributions, as shown in 
fig. ^: a sea contribution, where the sea quark or gluon is emitted from a soft Pomeron, and a 
valence contribution, where the valence quark is one of the three quarks of the nucleon. The precise 
microscopic structure of the soft Pomeron not being known, it is parameterized in the usual way 
a a Regge pole. 

Elementary nucleon-nucleon scattering can now be considered as a straightforward general- 
ization of photon-nucleon scattering: one has a hard parton-parton scattering in the middle, and 
parton evolutions in both directions towards the nucleons. We have a hard contribution Thard when 
the the first partons on both sides are valence quarks, a semi-hard contribution Tsemi when at least 
on one side there is a sea quark (being emitted from a soft Pomeron), and finally we have a soft 
contribution, when there is no hard scattering at all (see fig. We have a smooth transition from 




Figure 4: The soft elastic scattering amplitude T^oft (left), the hard elastic scattering amplitude Tjiard (middle) 
and one of the three contributions to the semi-hard elastic scattering amplitude Tgcmi (right). 

soft to hard physics: at low energies the soft contribution dominates, at high energies the hard 
and semi-hard ones, at intermediate energies (that is where experiments are performed presently) 
all contributions are important. 

Let us consider nucleus- nucleus {AB) scattering. In the Glauber-Gribov approach H, ^, the 
nucleus-nucleus scattering amplitude is defined by the sum of contributions of diagrams, corre- 
sponding to multiple elementary scattering processes between parton constituents of projectile 
and target nucleons. These elementary scatterings are exactly as discussed above, namely the 
sum of soft, semi-hard, and hard contributions: T2^2 — Tsoft + Tsomi + Tiiard- A corresponding 
relation holds for the inelastic amplitude T2^x- We use the above definition of a cut elementary 
diagram, which is graphically represented by a vertical dashed fine, whereas the elastic amplitude 
is represented by an unbroken line: 

= 72^2, ' = Z]x(^2^x)(72^x)*• 
This is very handy for treating the nuclear scattering model. We define the model via the elastic 
scattering amplitude Tab->ab which is assumed to consist of purely parallel elementary interactions 
between partonic constituents, as discussed above. The amplitude is therefore a sum of terms as 
the one shown in fig. ^. One has to be careful about energy conservation: all the partonic 
constituents (fines) leaving a nucleon (blob) have to share the momentum of the nucleon. So in 
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Figure 5: Elastic nucleus-nucleus scattering amplitude, being composed of purely parallel elementary interactions 
between partonic constituents of the nucleons (blobs). 

the explicit formula one has an integration over momentum fractions of the partons, taking care 
of momentum conservation. Having defined elastic scattering, inelastic scattering and particle 
production is practically given, if one employs a quantum mechanically self-consistent picture. Let 
us now consider inelastic scattering: one has of course the same parallel structure, just some of the 
elementary interactions may be inelastic, some elastic. The inelastic amplitude being a sum over 
many terms - Tab^x = "^ab^x ^ ^o be squared and summed over final states in order to 
get the inelastic cross section, which provides interference terms X]x(^AB^jf which 
can be conveniently expressed in terms of the cut and uncut elementary diagrams, as shown in fig. 
||. So we are doing nothing more than following basic rules of quantum mechanics. Of course a 
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Figure 6: Typical interference term contributing to the squared inelastic amplitude. 

diagram with 3 inelastic elementary interactions does not interfere with the one with 300, because 
the final states are different. So it makes sense to define classes K of interference terms (cut 
diagrams) contributing to the same final state, as all diagrams with a certain number of inelastic 
interactions and with fixed momentum fractions of the corresponding partonic constituents. One 
then sums over all terms within each class K , and obtains for the inelastic cross section 

aAB{s)^ I d26^r!M)(i^) 
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where we use the symbohc notation (fb = J cPbo J (P^Ba p{bA) J (P^bs pibs) which means inte- 
gration over impact parameter bo and in addition averaging over nuclear coordinates for projectile 
and target. The variable K is characterized by AB numbers mk representing the number of cut 
elementary diagrams for each possible pair of nucleons and all the momentum fractions and x~ 
of all these elementary interactions (so if is a partly discrete and partly continuous variable, and 
J2 is meant to represent J2 /)■ This is the really new and very important feature of our approach: 
we keep explicitly the dependence on the longitudinal momenta, assuring energy conservation at 
any level of our calculation. 

The calculation of ft actually very difficult and technical, but it can be done and we refer the 
interested reader to the literature 

The quantity fi^^'^^ (ii') can now be interpreted as the probability to produce a configuration 
K at given s and 6. So we have a solid basis for applying Monte Carlo techniques: one generates 
configurations K according to the probability distribution and one may then calculate mean 
values of observables by averaging Monte Carlo samples. The problem is that fl represents a 
very high dimensional probability distribution, and it is not obvious how to deal with it. We 
decided to develop powerful Markov chain techniques |Q| in order to avoid to introduce additional 
approximations. 

To summarize: We provide a new formulation of the multiple scattering mechanism in nucleus- 
nucleus scattering, where the basic guideline is theoretical consistency. We avoid in this way many 
of the problems encountered in present day models. We also introduce the necessary numerical 
techniques to apply the formalism in order to perform practical calculations. 

This work has been funded in part by the IN2P3/CNRS (PICS 580) and the Russian Foundation 
of Basic Researches (RFBR-98-02-22024). 
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